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The Multi-Band Template Analysis is a low-latency analysis pipeline for the detection of grav¬ 
itational waves to triggering electromagnetic follow up observations. Coincident observation 
of gravitational waves and an electromagnetic counterpart will allow us to develop a com¬ 
plete picture of energetic astronomical events. We give an outline of the MBTA pipeline, as 
well as the procedure for distributing gravitational wave candidate events to our astronomical 
partners. We give some details of the recent work that has been done to improve the MBTA 
pipeline and are now making preparations for the advanced detector era. 


1 Introduction 

Currently the LIGO ^ and Virgo ^ detectors are being brought back into operation after an 
extended period of upgrades and commissioning. This is an extremely exciting time in the 
gravitational wave (GW) community as we prepare for the beginning of the advanced detector 
era, when advanced LIGO ^ and advanced Virgo ^ come online. The advanced detectors will 
have a sensitive range for binary neutron stars (BNS) source a factor of 10 better than the 
initial detectors, which corresponds to an increase in the observable volume by three orders of 
magnitude. This will improve the “realistic” BNS detector rates from 0.02 yr“^ in the initial 
detector era to 40yr“^ in the advanced detector era^. 

The first scheduled observing run in the advanced detector era will be a three month run 
starting in September 2015 which will only include the two LIGO detectors while the Virgo 
detector finishes its upgrades and commissioning. Virgo will join the detector network for the 
first three detector observing run in 2016-2017 for six month. 

In this paper we will outline the key elements of the Multi-Band Template Analysis (MBTA) 
pipeline and the procedure for distributing GW candidate events for electromagnetic (EM) follow 
up observations, as shown in figure 1. In section 2.1 we explain how GW signals are extracted 
from the GW channel data of each detector. In section 2.2 we explain the criteria for selecting 
GW candidate events. In section 2.3 we give details of the post-MBTA event processing and 
explain how GW candidate events are distributed to our astronomical partners. In section 2.4 
we give details of some of the improvements recently implemented in MBTA. In sections we 
summarise the status of the MBTA pipeline in preparation for the advanced detector era. 

2 The Multi-Band Template Analysis 

MBTA^ is a low-latency coincidence analysis pipeline used to detect GWs from compact binary 
coalescences (CBCs). Some of the best understood and strongest sources of GWs for the LIGO 



4 - H1 
4-L1 
4- vi 
GW data 
generation 




Identify 
GW triggers 


Write to GW 
database 


Sky 

localisation 


Send alerts 
for EM 
follow up 


Figure 1 - Overview of the joint GW-EM observation pipeline. 


and Virgo detectors are CBCs consisting of two neutron stars (NS-NS), or a neutron star and a 
black hole (NS-BH). These systems have many possible mechanisms for producing EM counter 
parts The possibility of a strong GW signal and EM counterpart make these sources the focus 
of the MBTA pipeline, which aims to detect GW candidate events with low enough latency to 
trigger EM follow up observations by our astronomical partners^. 

2.1 Single detector analysis 

Each detector in the network is analysed independently, before the results are later combined 
to find GW candidate events. We obtain the calibrated GW channel data from each of the 
detectors, as well as basic data quality information, informing us of the status of the detectors. 
MBTA uses the standard matched filter^ to extract CBC signals from the GW channel data of 
each detector. To do this a bank of search templates is used to cover the parameter space of 
expected signals, this is generated at initialisation to keep the analysis latency as low as possible. 
This template bank is referred to as the “virtual” template bank, and covers the parameter space 
we are interested in. 

To reduce the computational cost of the matched filtering, which is the most computationally 
expensive element of the analysis, MBTA splits the matched filter across two (or more) frequency 
bands. The boundary frequency between the low and high frequency bands is selected so that 
the signal-to-noise ratio (SNR) is roughly equally shared between the two bands. This multi¬ 
band analysis procedure gives a reduction in the computational cost and we loose negligible 
SNR compared to a matched filter performed with a single band analysis. The reduction in 
computational cost comes from the fact that in each frequency band we can use shorter templates 
and so the phase of the signal is tracked over fewer cycles, this reduces the number of templates 
that is required to cover the equivalent mass space of a single band analysis. We also benefit 
from being able to use a reduced sampling rate for the low frequency band, which reduces the 
computational cost of the fast Fourier transforms involved in the filtering. 

Each frequency bands requires a separate “real” template bank, which is actually used to 
filter the data. For each template in the virtual template bank, a template from the low and 
high frequency real template banks are combined to reconstruct the result from a single band 
analysis. The real template banks and the parameters for combination and association with the 
virtual template bank are produced during initialisation to keep the analysis latency low. 

To further reduce the computational cost of the filtering, the template banks are split across 
multiple jobs and run in parallel across the parameter space. Once we have the results from 
each band across the full parameter space the results are coherently combined between the 
bands. Triggers are extracted from the match filter output in each band when SNR > 5, and a 
computationally inexpensive test is used to check the SNR distribution across the frequency 
bands is consistent with the expected signal SNR evolution. 

2.2 Coincidence events 

The trigger lists from the individual detectors are combined to find coincidence events. Time 
coincidence it checked using a simple time of flight consistency test between triggers in detector 
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Figure 2 - Skymap of GW candidate event G20190 


pairs. In the past we also used a mass coincidence criterion, but this has now been superseded by 
the exact match requirement. The exact match requires that triggers are found in all detectors 
with the same template parameters; the component masses and spins. The significance of each 
event is estimated by calculating the false alarm rate (FAR), the expected rate of coincidence 
triggers from noise only that have an equal or large SNR than the event. 


2.3 Event follow up 

GW candidate events found by the MBTA pipeline are sent to the Gravitational Wave Can¬ 
didate Event Database (GraCEDb) an automated archive for GW candidate events. When 
MBTA uploads a new event to GraCEDb rapid sky localisation is performed with Bayestar^^, a 
rapid Bayesian position reconstruction code, using the time, amplitude, and phase information 
reported by MBTA. Bayestar returns a probability skymap that is appended to the GraCEDb 
event and can be used to plan follow up observations around the most probably sky positions 
of the source. As well as sky localisation, any additional data quality vetoes can be applied to 
remove events that are associated with detector noise. 

In previous analysis runs^ GW candidate events have then undergone verification by a human 
monitor. The purpose of this human monitor was to review the events, consult the detector 
control rooms about the status of the detectors, and verify the data quality. In the past the entire 
pipeline, from data collection distributing GW candidate events to our astronomical partners, 
had a latency of 20-40 minutes; where the human monitor step gave the largest contribution to 
the latency. We give an example of a GW candidate event skymap^ in figure2. This event was 
distributed to our astronomical partners with a latency of 39 minutes and follow up observations 
were performed by Quest, ROTSE, SkyMapper, TAROT and Zadko. 


2.4 Recent improvements 

In preparation for the first observing run in the advanced detector era we have implemented a 
number of improvements to the MBTA pipeline. We are now able to run MBTA using spinning 
template banks, this was achieved by changing the interface for generating the template banks. 
Now that we are using the exact match, we can provide the coalescence phase of events for 
improved sky localisation with Bayestar. Finally, to remove noise events we have added a signal 
based consistency test. This uses the fact that a real CBC signal should produce a single loud 
peak in the matched filter output, whereas a noise trigger could possibly produce multiple loud 
peaks. Comparing the peak SNR of an event to the surrounding level we can veto events that 
do not behave like real CBC signals. 


3 Conclusion 


In this paper we have outlined the MBTA pipeline, and how the GW candidate events produced 
by MBTA are distributed to our astronomical partners for EM follow up observations. As was 
shown in previous observing runs, we are ready to perform this task and are now focusing on 
implementing improvements in preparation for the first advanced era observing run in September 
2015. We expect that as the detectors sensitivities are improved, we will soon be making the 
first GW detections and can begin to do real GW astronomy in coalition with our astronomical 
partners. 
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